Hadron-string cascade versus hydrodynamics in Cu + Cu collisions at ^/snn = 200 GeV 
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Single particle spectra as well as elliptic flow in Cu+Cu collisions at ^/sjvjv = 200 GeV are in- 
vestigated within a hadronic cascade model and an ideal hydrodynamic model. Pseudorapidity 
distribution and transverse momentum spectra for charged hadrons are surprisingly comparable 
between these two models. However, a large deviation is predicted for the elliptic flow. The forth- 
coming experimental data will clarify the transport and thermalization aspects of matter produced 
in Cu+Cu collisions. 



PACS numbers: 24.10.Lx,24.10.Nz,25.75.Ld,25.75.-q 

One of the primary current interests in the Relativis- 
tic Heavy Ion Collider (RHIC) experiments is to explore 
the properties of QCD matter far from stable nuclei, es- 
pecially the confirmation of the deconfmed and thcrmal- 
ized matter, i.e. the quark gluon plasma (QGP), which 
has been predicted from the lattice QCD calculations [1] . 
While high and medium p T observables such as the par- 
ton energy loss [2] and coalescence behavior of hadron 
elliptic flows [3] are generally believed to give strong evi- 
dences of high dense matter formation, hadrons at these 
momenta are not necessarily formed from thcrmalizcd 
matter. Therefore, low p T observables are also impor- 
tant to confirm whether equilibrium is achieved or not. 

Elliptic flow [4] is one of the promising observables 
to study the degree of thermalization for QCD matter 
produced in heavy ion collisions since it is believed to 
be sensitive to the properties of the matter at initial 
stages and the collision geometry [5, 6]. Indeed, inci- 
dent energy as well as impact parameter dependences 
of elliptic flow have been investigated extensively. El- 
liptic flow, i.e. the momentum anisotropy with respect 
to the reaction plane V2 = (cos(20)), has been mea- 
sured in a wide energy range from GSI-SIS (E[ nc < 1A 
GeV) [7], BNL-AGS (E inc = 2-114 GeV) [8], to CERN- 
SPS (.Bine = 40 - 158A GeV) [9], in addition to BNL- 
RHIC [10]. Measured collective flows are well reproduced 
by nuclear transport models assuming the momentum de- 
pendent nuclear mean- field at SIS to AGS (-Einc — 0.2 — 
ITAGeV) [11, 12] and SPS (E inc = 40, 158AGcV) [13] 
energies, whereas elliptic flow at RHIC at mid-rapidity 
is underestimated in noncquilibrium transport models 
which do not include explicit partonic interactions [14- 
16]. It is also reported that hadronic models explain 
elliptic flow only at low transverse momentum pt < 1 
GeV/c at RHIC [17]. Partonic interactions followed by 
quark coalescence hadronization mechanism are proposed 
in Rcf. [18] to account for the experimental data on el- 
liptic flow. Note, however, that hadronic cascade models 
reproduce elliptic flow in forward/backward rapidity re- 
gions at RHIC [19]. 

On the other hand, in Au+Au collisions at RHIC ener- 



gies the magnitude of V2 and its transverse momentum pt 
and mass m dependences are close to predictions based 
on ideal and non-dissipativc hydrodynamics simulations 
around midrapidity (| 7/ | ;$ 1), in the low transverse mo- 
mentum region (pr~l GeV/c), and up to semicentral col- 
lisions 5 fm) [20, 21]. This is one of the main results 
which leads to a recent announcement of the discovery 
of perfect fluidity at RHIC [22]. (See Ref. [23] for recent 
rcintcrpretation of the RHIC data based on current hy- 
drodynamic results.) Despite the apparent success near 
midrapidity at RHIC, ideal hydrodynamics overestimates 
the data at lower incident energies (SIS, AGS and SPS) 
as well as in forward/backward rapidity regions at RHIC 
probably due to the lack of dissipative effects. 

We study Cu+Cu collisions at RHIC in the present 
work, which is a complementary study of elliptic flow 
in Au+Au collisions. The particle density and the size 
of the system are smaller in Cu+Cu collisions than in 
Au+Au collisions. So the reasonable agreement of hy- 
drodynamic results with Au+Au data may be spoiled 
in Cu+Cu collisions and a non-equilibrium hadronic de- 
scription can be relatively important even at RHIC en- 
ergies. Therefore we employ both a hadronic transport 
model JAM and a hydrodynamic model to make pre- 
dictions for elliptic flow in Cu+Cu collisions. Below we 
briefly summarize hadron-string cascade JAM [24] and a 
hydrodynamic model [25] adopted in this paper. 

A hadronic transport model JAM simulates nuclear 
collisions by the individual hadron-hadron collisions. 
Soft hadron productions in hadron-hadron scattering are 
modeled by the resonance and color string excitations. 
Hard partonic scattering is also included in line with HI- 
JING [26]. Color strings decay into hadrons after their 
formation time (r ~ 1 fm/c) according to the Lund string 
model PYTHIA [27]. Hadrons within their formation 
time can scatter with other hadrons assuming the ad- 
ditive quark cross section. This simulates constituent 
quark collisions effectively which is known to be impor- 
tant at SPS energies [28]. Therefore, matter initially cre- 
ated in collisions is represented by the many strings at 
RHIC, which means that there is no QGP in the model. 
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Default parameters in JAM are adopted in this work 
except for a little wider p T width in the string decay 
and a larger partonic minimum p T (po = 2.7GeV/c) 
to fit charged hadron p T spectrum in pp collisions at 
y/sNN = 200 GeV. In addition to hadron-hadron col- 
lisions, nuclear mean field is incorporated in JAM and 
its effects are known to be important at AGS and SPS 
energies [13], but mean field is not expected to play ma- 
jor roles at RHIC. We have thus neglected nuclear mean 
field in this work. The detailed description of JAM can 
be found in Ref. [24]. 

Two of the authors (T.H. and Y.N.) have already de- 
veloped another dynamical framework to describe three 
important aspects of relativistic heavy ion collisions [25] , 
namely color glass condensate (CGC) for collisions of two 
nuclei [29-31], hydrodynamics for space-time evolution 
of thermalized matter [32], and jet quenching for high 
Pt non-thermalized partons [33] . Along the line of these 
works, we use the same model in this study. However, our 
aim is to study the bulk properties of matter produced in 
Cu+Cu collisions. In this paper, we neither include jet 
components in this model nor discuss jet quenching, un- 
like a series of the previous work [34] . So hydrodynamic 
results to be presented below include purely boosted ther- 
mal components without any semi-hard components. 

In Ref. [25], a systematic hydrodynamic analysis in 
Au+Au collisions at ^snn = 200 GeV was performed by 
using initial conditions taken from the CGC picture for 
the colliding nuclei. In the conventional hydrodynamic 
calculations, one chooses initial condition for hydrody- 
namic equations and thermal freezeout temperature T th 
so as to reproduce the observed particle spectra, such as 
(pseudo)rapidity distribution and transverse momentum 
distribution. So it is believed that hydrodynamics has 
a less predictive power compared with cascade models. 
However, if the initial particle production at high colli- 
sional energies is supposed to be universal as described 
by the CGC, hydrodynamics with CGC initial conditions 
can predict particle spectra. Here, we employ the IC-n, 
i.e. a prescription that the number density produced in 
a CGC collision is matched to the hydrodynamic initial 
condition [25], to obtain the initial distribution of ther- 
modynamic variables at the initial time to- Once the 
initial condition is obtained, one solves hydrodynamic 
equation d^T^ = in the three-dimensional Bjorken 
coordinate (r,r] s ,x,y) [21]. Here we neglect a dissipativc 
effect and a finite (but probably tiny) baryon density. As- 
suming N c = Nj = 3 massless partonic gas, an ideal gas 
equation of state (EOS) with a bag constant B 1 / 4 = 247 
McV is employed in the QGP phase (T > T c = 170 
MeV). We use a hadronic resonance gas model with all 
hadrons up to A(1232) mass for later stages (T < T c ) 
of collisions. We take into account chemical freezeout 
separated from thermal freezeout [35] as required to ob- 
tain sufficient yields for heavier particles. Specifically, we 
assume that chemical freezeout temperature T ch = 170 
MeV and kinetic freezeout temperature T th = 100 MeV. 
Note that the slope of px spectra becomes insensitive to 
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FIG. 1: (Color online) Pseudorapidity distributions for 
charged hadrons in Cu+Cu collisions at ^/snn = 200 GeV 
for impact parameters 6=1,2, and 5 fm. Circles correspond 
to the result of JAM. Lines denote the results of hydrody- 
namics for T th = 100 MeV. 



T th (while V2 becomes sensitive) when chemical freezeout 
is taken into account [35]. In the calculation of V2{r]), we 
also use T th = 160 MeV for comparison. If the strongly 
coupled QGP (sQGP) core expands as a perfect fluid and 
the hadronic corona does as a highly dissipativc gas as 
suggested in Ref. [23], the resultant ^(f?) and V2(pr) are 
expected to be frozen after hadronization [36, 37] due to 
the strong viscous effect. Moreover, T th should be higher 
for a smaller size of the system [38] as observed in the cen- 
trality dependence of the pr spectra [39] . So the freeze- 
out picture in Cu+Cu collisions can be different from that 
in Au+Au collisions. In the following predictions for hy- 
drodynamic elliptic flow, we show the results for T th = 
100 and 160 MeV. For further details of the hydrody- 
namic model used in this work, see Refs. [25, 34, 35]. 

We first compare the bulk single particle spectra be- 
tween JAM and hydrodynamics. We emphasize again 
that our hydrodynamic results are insensitive to a choice 
of T th for transverse and rapidity distributions of charged 
hadrons. We show results of the pseudorapidity distribu- 
tion dN/drj for charged hadrons in Fig. 1 at impact pa- 
rameters 6=1,2 and 5 fm. It is seen from this figure that 
the shape and the magnitude of the distributions from 
JAM are almost similar to those from hydrodynamics. 

In Fig. 2, we compare JAM and hydrodynamic results 
of the p T spectra for charged hadrons at impact param- 
eters of 6 = 1,2 and 5 fm for |t;| < 0.33. Accidentally, 
these results agree well with each other in transverse mo- 
mentum range of p T < 2 GcV/c. Deviation at higher 
transverse momentum is due to the lack of jet compo- 
nents in the hydrodynamic simulations. 

At least within our models, two distinct pictures, i.e. 
pictures of coherent particle production via CGC com- 
bined with sequential sQGP expansion and of trans- 
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FIG. 2: (Color online) Comparison of the transverse momen- 
tum distributions for charged hadrons between JAM (circles) 
and hydrodynamics for T th = 100 MeV (lines) at \r)\ < 0.33 
in Cu+Cu collisions at ^/snn = 200GeV. The results for 6=2 
fm (5 fm) are scaled by 10 -2 (10 -4 ). 




FIG. 4: (Color online) Elliptic flow V2 for charged hadrons as a 
function of pseudorapidity r\ in Cu+Cu collisions at ^/snn = 
200 GeV at impact parameter b = 5 fm. Circles connected 
by line show the results of JAM. Triangles and squares with 
lines show the results of hydrodynamics with T th =100 MeV 
and 160 MeV, respectively. 




01 2345678 
impact parameter b(fm) 



FIG. 3: (Color online) Elliptic flow 112 for charged hadrons 
at mid-rapidity as a function of impact parameter in Cu+Cu 
collisions at ^/snn = 200GeV. Circles connected by line show 
the results of JAM. Triangles and squares with lines show the 
results of hydrodynamics with T th =100 MeV and 160 MeV, 
respectively. 



ports of secondary hadrons after hadron-hadron colli- 
sions summed up by an overlap region of colliding nuclei, 
are indistinguishable in the bulk single hadron distribu- 
tions in Cu+Cu collisions. Note that free parameters in 
the "CGC+hydro" model has been fixed by fitting the 
charged multiplicity in Au+Au collisions at midrapidity. 
We also note that parameters in JAM are already fixed 
to fit the data in pp collisions at = 200 GeV. 

In Fig. 3, we show the impact parameter b dependence 
of the elliptic flow V2 at mid-rapidity for charged hadrons. 



In the hydrodynamic calculations, kinetic freezeout tem- 
peratures T th = 100 MeV and 160 MeV are chosen. 
While single particle spectra from JAM and hydrody- 
namics look very similar, a clear difference of V2(b) is 
seen: v-i grows almost linearly with b in hydrodynamics, 
which is the same as the case in Au+Au collisions, while 
we find a peak at around b = 6 fm in JAM and that the 
magnitude is only around 20% of the hydrodynamic pre- 
diction with T th = 100 MeV. The two distinct pictures 
within our approach appear differently in the ccntrality 
dependence of elliptic flow. Due to the smaller initial 
energy density in Cu+Cu collisions compared to Au+Au 
collisions, the spatial anisotropy is still out-of-plane just 
after the hadronization and V2 continues to be generated 
even in the late non-viscous hadronic stage in the ideal 
hydrodynamic simulation. The data is expected to be 
comparable with the result for T th = 160 MeV if the 
initial energy density is large (erj 3> 1 GcV/fm 3 ) and 
the equilibration time is small (to ~ 1 fm/c) enough to 
create the sQGP phase in Cu+Cu collisions. On the 
contrary, one expects that it takes more time to reach 
equilibrium (to > 1 fm/c) and that the system may not 
reach the equilibrated sQGP state since the system size 
and the produced particle number are small compared 
with those in Au+Au collisions. In that case, the data 
will be comparable with the result from JAM. 

Pseudorapidity dependences of the elliptic flow from 
JAM and hydrodynamics are compared with each other 
to understand the longitudinal dynamics in Cu+Cu col- 
lisions in Fig. 4. In JAM, we find almost flat behavior 
of 1^2(17 ) around midrapidity (\r)\ < 2), where the charged 
hadron i] distribution also shows flat behavior. In JAM, 
elliptic flow is slowly generated (t < lOfm) as the hadrons 
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FIG. 5: (Color online) The calculated elliptic flow Vi of 
charged hadrons as a function of transverse momentum pr 
for Cu+Cu at y^sjvjv = 200 GeV for different impact param- 
eters, b = 1 fm (top), 2 fm (middle), and 5 fm (bottom). 
Circles connected by line show the results of JAM. Triangles 
and squares with lines show the results of hydrodynamics with 
T th =100 MeV and 160 MeV, respectively. 



are formed from strings after some formation times. In 
the hydrodynamic calculations, we show the results for 
T th = 100 and 160 MeV, which could be an upper and 
a lower limit of the ideal hydrodynamic prediction re- 
spectively. V2{rj) for T th = 100 MeV becomes a trape- 
zoidal shape, which looks similar to the result in the pre- 
vious hydrodynamic study in Au+Au collisions [21, 35]. 
V2{rj) from ideal hydrodynamics for T th = 160 MeV is 
also shown as a possible result for the situation [23] in 
which V2 is generated by the perfect fluid of the sQGP 



core and is not generated significantly in the dissipative 
hadronic corona like the result from JAM. Indeed, V2(r]) 
for T th = 160 MeV appears to be a triangle shape which 
looks similar to the shape in Au+Au data observed by 
PHOBOS [40]. 

In Fig. 5, we compare transverse momentum p T depen- 
dence of elliptic flow for charged hadrons. Hydrodynamic 
predictions are of course larger than the ones of JAM. In 
JAM, V2 starts to be saturated at around 0.8 GeV, and 
the behavior is qualitatively similar to that in Au+Au 
collisions and another theoretical prediction in Cu+Cu 
collisions [41]. It should be noted that we will also find a 
mass dependent saturating behavior of t>2 (pt) when semi- 
hard components are combined with the hydrodynamic 
components [25]. 

In summary, we have investigated low-p T observables 
in a hadron-string cascade model JAM [24] and a hydro- 
dynamical model [25] in Cu+Cu collisions at ^/s^n = 
200 GeV. For dN/dij and p T -spectra for charged hadrons, 
we have obtained good agreement between JAM and hy- 
drodynamics. However, clear deviations between model 
predictions are found in V2 as a function of central- 
ity, pseudorapidity 77, and transverse momentum p T for 
charged hadrons. The lack of elliptic flow in hadronic 
transport models compared to the ideal hydrodynamic 
predictions is due to the initial particle production being 
performed by string decays which only generate a limited 
amount of transverse momentum for the produced par- 
ticles in conjunction with the formation time for these 
hadrons. In that sense, the "EOS" of hadron-string cas- 
cade models in the very early stage is to be considered as 
a super-soft one and cannot generate sufficient pressure 
needed for elliptic flow to develop. In addition, even if full 
thermalization is achieved in the hadron cascade model, 
a higher viscosity in the hadron cascade model would 
yield a lower elliptic flow than the ideal non-viscous hy- 
drodynamics in the hadron phase. Therefore, in order to 
interpret correctly the result, V2, cascade < V2, hydro, which 
has been seen also in SPS energies, we should study the 
dissipative effects carefully Measurements of pseudora- 
pidity and transverse momentum dependence of elliptic 
flow (v2(b), 1*2(77), and ^(pt)) m Cu+Cu collisions at 
RHIC will provide very important information for trans- 
port aspects of QCD matter in heavy ion collisions. 
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